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Abstract 

Recently, a new mechanism, which explains why the three gauge coupling con- 
stants meet at a certain scale in the minimal super symmetric standard model, has 
been proposed in a scenario of grand unified theories with anomalous U(1)a gauge 
symmetry. It is a non-trivial result that although there are many superheavy fields 
whose mass scales are below the unification scale, this mechanism explains this fact 
by means of one-loop renormalization group equations. Since the unification scale 
generically becomes below the usual GUT scale, 2 x 10 16 GeV, and proton decay 
via dimension 5 operators is suppressed, the scenario predicts that proton decay via 
dimension 6 operators, p — ► e7r°, will be observed in the near future. In this paper, 
we attempt to estimate a reasonable range of values of the lifetime of the proton 
predicted within this scenario by using a two-loop renormalization group calculation 
and the ambiguities of 0(1) coefficients. 
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1 Introduction 



Low energy supersymmetry (SUSY), which was originally introduced for stabilization of 
the weak scale, plays a critical role in explaining the hierarchical gauge couplings of the 
standard model in the context of SU(5) grand unified theory (GUT). It is a non-trivial 
result that the three gauge couplings meet at a GUT scale ~ 2 x 10 16 GeV if a rea- 
sonable SUSY breaking scale is assumed. However, the experimental lower bound of the 
proton lifetime via dimension 5 operators has become so severe that many GUTs that 
naturally realize coupling unification, for example, the minimal SU(5) GUT, have been 
rejected. [[I], |2|, [| This is a kind of puzzle in the SUSY GUT scenario. Of course, if we 
consider the coincidence of the three gauge couplings at some scale to be accidental, then 
this puzzle is not difficult to solve. For example, after suppression of dimension 5 proton 
decay, the gauge couplings can be caused to meet at a certain scale by tuning the mass 
scales under the GUT scale. However, if the coincidence of these couplings is not acci- 
dental, solving this puzzle is not so simple, although some solutions have been proposed. 
Most of these solutions aim to show that the minimal SUSY standard model (MSSM) is 
realized under the GUT scale with the dimension 5 proton decay suppressed or for- 
bidden. There are several solutions of this type in the context of extra dimensions. One 
employs parity assignment, which forbids dimension 5 operators. || [jj The other employs 
wave function suppression due to localization of quark and lepton fields. [§, || Even in 
the context of four-dimensional field theory, solving this puzzle is possible by introducing 
a special vacuum structure of one or two adjoint Higgs fields of SO (10). |T(| However, 
recently, another type of solution has been proposed JXTJ] in the context of GUT with 



anomalous U(l) gauge symmetry, ||12|| whose anomaly is cancelled by the Green-Schwarz 
mechanism.^] It is surprising that in the GUT scenario with a simple unification group, 
gauge coupling unification is realized generically, despite the fact that many superheavy 
fields, whose mass spectrum does not respect the GUT symmetry, become lighter than 
the unification scale A^ and there are several gauge symmetry breaking scales. This is 
because the mass spectrum of superheavy fields and the symmetry breaking scales are 
determined by the anomalous U(1)a charges, and most of the charges are cancelled un- 
der the conditions of coupling unification. (The unique exception is the charge of the 
doublet Higgs.) Moreover, the GUT scenario has many interesting features which we 



now describe. [p| [14|, 15, [16], [17]] 1) The interaction is generic in the sense that all the 



interactions that are allowed by the symmetry are introduced. Therefore, once we fix the 
field content with their quantum numbers (integers), all the interactions are determined, 
except the coefficients of order 1. 2) It naturally solves the so-called doublet-triplet (DT) 



splitting problem, [IS] using the Dimopoulos-Wilczek (DW) mechanism. [19, 3) It re- 



produces the realistic structure of the quark and lepton mass matrices, including neutrino 



bi-large mixing, pl[] using the Froggatt-Nielsen (FN) mechanism. p2[ 4) The anomalous 
U(1)a accounts for the hierarchical structure of the symmetry breaking scales and the 
masses of heavy particles. 5) All the fields, except those of the minimal SUSY standard 
model (MSSM), can become heavy. 6) The gauge couplings are unified just below the 



1 Recently, it was pointed out that allowing arbitrary soft masses and fermion and sfermion mixing, 
there is still only a small portion of parameter space consistent with experimental results in the 'decou- 
pling' region. H 
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usual GUT scale, A G ~ 2 x 10 16 GeV. 7) In spite of the lower unification scale, proton 
decay via dimension 6 operators, p — > e + n°, is still within the experimental bound and 
therefore we expect to observe proton decay in the near future. 8) The cutoff scale is 
lower than the Planck scale. 9) The /i problem is also solved. 

In the above-mentioned scenario, one of the most interesting predictions regards pro- 
ton decay. Because the dimension 5 operators are suppressed, the main decay mode of 
proton decay is due to dimension 6 operators. Therefore, a more accurate estimate of 
the unification scale or the cutoff scale is important in order to obtain a more accurate 
prediction of the lifetime of the proton. In this paper, we determine an allowed range of 
values of the cutoff scale for several GUT models by means of two-loop renormalization 
group equations (RGEs) and using the freedom of 0(1) coefficients. 

2 Gauge coupling unification (one-loop renormal- 
ization group) 

In our scenario, because generic interactions are introduced, the order of every coefficient is 
determined by anomalous U(\)a charges. Therefore, the gauge symmetry breaking scales 
and the mass spectrum of superheavy fields are also determined by these charges. Thus 
we can examine whether the gauge couplings meet at the GUT scale once we determine 
all the charges. This is a consistency check of our scenario, and it has been shown that it 
is realized in a non-trivial way. In this section, we review this point using one-loop RGEs. 

First, we note that the symmetry breaking scales are determined by anomalous U{1)a 
charges. Generically, the vacuum expectation value (VEV) of the gauge singlet operator 
O is determined by its charge o as 

(O) ~ A"°. (1) 

(Here A is the ratio of the cutoff scale A to the VEV of the Froggatt-Nielsen field 0. In 
our scenario, we adopt A ~ 0.22. ) Actually, in our scenario, the VEV of the adjoint field 
A of 5*0(10) becomes of the Dimopoulos-Wilczek type : 

(A) = ir 2 x diag(-u, v, v, 0, 0). (2) 

This fact plays an important role in realizing doublet-triplet splitting, and the scale is 
determined by the charges a as v ~ A~ a A. Throughout this paper we denote all superfields 
by uppercase letters and their anomalous U{1)a charges by the corresponding lowercase 
letters. We often use units in which A = 1. The VEV of A breaks 5O(10) into SU(3) C x 
SU(2) L x SU(2) R x U(1) B -l- The VEVs of the spinors C(16) and £(16), which break 
SU{2) R x U(1) B -l into U(1) Y , are obtained from (CC) ~ A~( c+5 ) and the Matness 
condition | (C) \ = \ (C"j | as 

| (C) | = | (C) | ~ A^ c+C ~)/ 2 . (3) 

Because | (C) \ < \ (A) \ is required to realize doublet-triplet splitting, at the scale A^ = 
(A) ~ A~ a , the 5O(10) gauge group is broken into SU(3) C x SU(2) L x SU(2) R x U(1) B -l, 
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which is broken into the standard gauge group SU(3)c x SU(2) L x U(1)y at the scale 
A c = (C) ~ A~( c+C ")/ 2 . 

Second, we explain how the mass spectrum of superheavy fields is determined by 
anomalous U(1)a charges. Using the definitions of the fields Q(3, 2) 1 , U c (3, l)_a , -D c (3, 1) 1 

6 3 3 

L(l,2)_i, i? c (l, l)i, iV c (l, l)o and Y(3, 2)_b, along with their conjugate fields, and 

2 6 

G(8, l)o and W(l,3)o with the standard gauge symmetry, under SO (10) D SU(5) D 
SU(3) C x 5C/(2) L x U(1) Y , the spinor 16, vector 10 and adjoint 45 of SO(10) are de- 
composed as 

16 -> [Q + f/ c + £ c ] + [D c + L] + JNT, (4) 

10 5 1 

10 -> [D C + L] + [D C + Z], (5) 
. ' „ ' 

5 5 

45 -> [G + IU + X + X + N c ] + [Q + U c + E c ] + [Q + U c + E c ] + N c . (6) 

S v ' S v ' S v ' 

24 10 TO 

A straightforward calculation of the mass matrices Mj of the superheavy fields I = Q, 
U c , E c , D c , L, G, W and X shows that 

detM/ = A^ c % (7) 

where the quantities q are the anomalous U(1)a charges of the superheavy fields. 

Third, we carry out an analysis based on the RGEs up to one loop. The conditions of 
gauge coupling unification are 

5 

a 3 (A A ) = a 2 (A A ) = ^a Y (A A ) = «i(Aa), (8) 

2 

where af 1 ^ > A c ) = |o;^ 1 (/i > Ac) + f^s-zXA 4 > ^-c)- Here ax = f^, and the pa- 
rameters gx (X = 3, 2, R, B — L, Y) are the gauge couplings of SU (3)c, SU (2) l, SU (2)r, 
U(1)b-l and U(l)y, respectively. 

The gauge couplings at the scale Aa are roughly given by 

1 (, , (M SB \ ^— v / vrii \ 12, (A c 



a-(A A ) = «;\m sb) + - Mn^j + EA^m^j--m^ 



(9) 

a 2 \A A ) = a 2 - 1 (M 5B ) + i-[6 2 ln(^)+EA6 2 an(^)), (10) 



a^(A A ) = «3 1 (^) + ^(^ln(^)+EA^ln(^)j, (11) 

where Msb is a SUSY breaking scale, (61,62,^3) — (33/5, 1, —3) are the renormalization 
group coefficients for the minimal SUSY standard model (MSSM), and Ab a i (a = 1,2,3) 
are the corrections to the coefficients from the massive fields with mass m,-. The last term 
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in Eq. (g) is from the breaking SU(2) R x U(1) B -l -> U(1) Y caused by the VEV (C). 
Because the gauge couplings at the SUSY breaking scale M SB are given by 



a T \M SB ) = a- G \K G ) + L{, M {^ 



(i = 1,2,3) 



(12) 



where a G 1 (Ac) ~ 25 and A G ~ 2 x 10 16 GeV, the above conditions for unification can be 
rewritten as 



,'7 



,'7 



,'7 



5 U G 



(13) 
(14) 
(15) 



Here, 77 represents the ranks of the mass matrices of the superheavy fields, Mj. 
corrections to the renormalization coefficients Ab a j are given in the following table: 
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Q + Q 


U c + U c 


E c + E c 


D c + D c 


L + L 


G 




x + x 
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8 

R 


6 
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2 

R 


3 

R 








5 


Ab 2I 


3 











1 





2 


3 


Ab 3I 


2 


1 





1 





3 
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The 



The unification conditions cti(A^) = a 2 (A y i), ai(A^) = 0:3 (A^) and a 2 (A J 4) = o 3 (Aa) 
be rewritten as 



can 



A G \ D ( detM L \ /detM Q \ 4 / det Mq \ 3 / det M w \ 5 



A G 7 VA A 7 VdetM DCy / ^detMc/y/ ^detA^,/ ^detM^ 

.-f D c+f L -4f u c-3f E c+7fQ-5f x +5f w 
— ly A i 

Ka \ W f Ac \ 4 / det Mjyc \ ( det M Q \ / det M \ 2 /det M G \ 



A 



G 



A. 



det M L I I det Afy / V det M £c / I det M x 



A 



A^ \ 4 / det Mjyc \ fdetMu\ /detM G \ /detM G x 



AqJ \ det M L J \det M Q J \det M w J \ det M 



x 



A 



-f£ +r D c-r Q +ru~2r w -r x +3r G 



(16) 



(17) 



'1? 



Note that the above conditions are dependent only on the ratio of the determinants of 
the mass matrices that are included in the same multiplet of SU(5) and on the symmetry 
breaking scales A^ and A G . If all the component fields in a multiplet were superheavy, 
the above ratios would be of order 1. However, because some of the component fields, for 
example massless Higgs doublets or Nambu-Goldstone modes, do not appear in the mass 
matrices generically, the above ratios are dependent only on the charges of these massless 



4 



modes. If all the fields other than those in MSSM become superheavy, the above ratios 
are easily estimated as 

, A -(ft«+h d ) j ( 19 ) 

> A c+5 " 2a , (20) 

> X~ 2a , (21) 





detM L 




det M D c 


detMg 


det M uc 


det M E - 


det M E c 


detM G 


detM w 


detMx 


detMx 



where h u and are the anomalous U(1)a charges of the massless Higgs doublets H u and 
Hdi respectively. Then the conditions for coupling unification become 

A ~ A^A G , (22) 
A ~ A — —A G , (23) 
A ~ A r^A G . (24) 

Therefore the unification conditions become h u + ha ~ 0, and thus the cutoff scale must 
be taken as A ~ A G . It is obvious that if the cutoff scale were some other scale (for 
example, the Planck scale), then in MSSM the three gauge couplings would meet at that 
scale. This implies that in this scenario it is not accidental that the three gauge couplings 
meet at some scale in MSSM. Note that these results are independent of the details of 
the Higgs sector, because the requirement that all the fields other than those in MSSM 
become superheavy determines the field content of the massless fields, whose charges are 
important to examine whether the gauge couplings meet at the unification scale A^. The 
above argument can also be applied to the scenario of Eq unification, though instead of 
the usual doublet Higgs charge h, we have to use effective Higgs charges, 

h ef[ = h+^((f ) -4>), (25) 



where E 6 is broken into 50(10) by the VEV | ($) | = | ($) | ~ A"3(*+*). 

Note that the condition h ~ does not mean that h = 0, because there is an ambiguity 
involving order 1 coefficients, and we have used only one-loop RGEs. However, the above 
analysis is useful in providing a rough picture of the behavior in which we are interested. 

3 Proton decay 

In order to understand how proton decay via dimension 5 operators is suppressed in our 
scenario, we now examine the mass matrix of triplet Higgs. From the interaction 

W = \ h+h ' +a HAH' + \ 2h 'H' 2 , (1) 

the mass matrices of doublet and triplet Higgs are 

(5 H ,5 H ,)l xh+h , +a{A) x2h , l(= I. (2) 
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Here H and H' represent 10 of 50(10) and their charges h < and h! > 0. Note that 
the H 2 term is forbidden by the SUSY zero mechanism. The colored Higgs obtain masses 
of order \ h+h +a (A) ~ \ h + h Because in general, \ h+h > \ 2h , proton decay is naturally 
suppressed. The effective colored Higgs mass is estimated as mf^ ~ (X h+h ) 2 /X 2h = X 2h , 
which is larger than the cutoff scale, because h < 0. If the cutoff scale is the usual GUT 
scale, A G = 2x 10 16 GeV, the condition mf > 10 18 GeV requires h < — 2. 

Strictly speaking, the condition for gauge coupling unification, h = 0, is not satisfied. 
However, we emphasize that the ambiguity of the coefficients of order 1 can naturally allow 
for the recovery of coupling unification. In order to change the unification condition, we 
have to consider the terms including the adjoint field A whose VEV (A) breaks GUT 
symmetry. In addition to the mass term m^II, higher- dimensional operators ^XAX 
must be taken into account. In the usual GUT scenario, such a correction is much smaller 
than the tree level mass term, because of the suppression factor Therefore, even if 
m ~ rrix < (A) is realized due to some symmetry, it is not easy to change the unification 
condition without finetuning or some special assumptions, for example, forbidding the 
mass term by some symmetry. However, the GUT with anomalous U(1)a symmetry 
naturally realizes mx ~ m \i because mx ~ X 2x , m ~ \ 2x + a anc l (A) ~ A _a . Moreover, 
in addition to these terms, XA n X also contributes to the mass of the X field after 
developing the VEV (A) if these terms are allowed by the symmetry. In other words, 
almost none of the coefficients respect GUT symmetry. Therefore, by adjusting such 
coefficients we can naturally change the unification condition. 

Finally, we recall proton decay via dimension 6 operators in our scenario. Because the 
cutoff scale A is around the usual GUT scale, = 2 x 10 16 GeV, and the unification 
scale becomes A _a (a < 0), proton decay via dimension 6 operators may be seen in 
future experiments. If we roughly estimate the lifetime of the proton using the formula in 
Ref. Q and a recent result provided by a lattice calculation for the hadron matrix element 
parameter we obtain 

This value! is near the present experimental limit. p3 



4 Two-loop analysis 

As mentioned above, in our scenario, the cutoff scale A is deeply related to the 5O(10) 
breaking scale (A) ~ A _a , and also to the mass of the massive gauge bosons, which cause 
proton decay. Therefore, to estimate the lifetime of the proton, we need to know how 
large the cutoff scale is. In addition, to realize somewhat accurate predictions, a two- 
loop analysis is necessary, because in our scenario there are many rather light superheavy 
particles, and therefore the coupling constants become larger than those of MSSM. In 

2 Note that this value is independent of the gauge coupling 510 at the unification scale. This is 
because the ratio ~§- is independent of gio , because the mass of massive gauge bosons mx is given by 

m x = \/2giov. 
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Table 1: Anomalous U(1)a charges of the non-trivial representational Higgs fields 
(A(45, -), A'(45, -), C(16, +), C(16, +), C'(16, -), C'(16, -), #(10, +) and #'(10, -)) 
and of the matter fields (^i(16, +), tf 2 (16, +), * 3 (16, +) and T(10, +)) in some models 
of the SO(10) scenario. Here, the sign + or — represents the Z 2 parity, and we assign 
odd R-parity for the matter fields. 





A 


A' 


C 


C 


a 


c 


H 


H' 


#1 


^2 


^3 


T 


i 


-1 


3 


-4 


-1 


3 


6 


-3 


4 


9/2 


7/2 


3/2 


5/2 


ii 


-1 


3 


-3 





2 


5 


-2 


3 


4 


3 


1 


2 


iii 


-1 


3 


-4 


-1 


3 


6 


-4 


5 


5 


4 


2 


3 


iv 


-1 


3 


-7/2 


1/2 


3/2 


11/2 


-3 


4 


9/2 


7/2 


3/2 


2 


V 


-1 


3 


-1 


-2 


4 


3 


-6 


7 


6 


5 


3 


6 


vi 


-1/2 


3/2 


-4 


-1 


5/2 


11/2 


-3 


7/2 


9/2 


7/2 


3/2 


5/2 


vii 


-1/2 


3/2 


-1 


-2 


7/2 


5/2 


-6 


13/2 


6 


5 


3 


6 



Table 2: Anomalous charges of the non-trivial representational Higgs fields 

(A(78, -), A'(78, -), $(27, +), $(27, +), C(27, +), (7(27, +), C"(27, -) and C"(27, -)) 
and of the matter fields (*i(27, +), * 2 (27, +) and ^ 3 (27, +)) in some models of the E 6 
scenario. Here, the sign + or — represents the Z 2 parity, and we assign odd R-parity for 
the matter fields. The MSSM doublet Higgs are contained in $. 





A 


A' 


$ 


$ 


C 


C 


O 


C" 


*i 


^2 


*3 


I 


-1/2 


5/2 


-3 


2 


-5 


-1 


13/2 


13/2 


9/2 


7/2 


3/2 


11 


-1/2 


5/2 


-3 


1 


-4 


-1 


13/2 


11/2 


9/2 


7/2 


3/2 


m 


-1 


4 


-3 


2 


-6 


-2 


7 


8 


9/2 


7/2 


3/2 



fact, in the E§ scenario, they sometimes become too strong to rely on the perturbation 
analysis. However, in the SO (10) scenario, the two-loop effect is still not excessively large, 
and we can treat it as a small correction to the one-loop analysis. 

In this section, for several models, we investigate how large the cutoff scale can be by 
using the two-loop RGEs and altering 0(1) coefficients in the range of < y < 7/ max 
as free parameters. The detailed procedure of this analysis is explained in Appendix [A|. 

4.1 Models 

There are several models that may work well. Among them we examine some represen- 
tative models in Tables [1] and 0. 

4.1.1 50(10) scenario 

There are four groups of models, characterized by different values of h: ii with h = —2; i,iv 
and vi with h = —3; iii with h = —4; v and vii with h = —6. A model with a larger charge 
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Table 3: A max in units of 10 16 GeV for some models of the SO(10) scenario. Basically, 
the parameters are set as y max = 2, A = 0.22 and a~ 1 (Mz) = 8.44, but in each case, one 
of these values is different from these, as indicated in the first row. The notation " — " 
means that there is no solution for gauge coupling unification. 







v - 2 2 


2/ max — 2 2 


A = 0.20 


A = 0.25 


a~ l = 8.30 


a; 1 = 8.58 


1 loop 


i 


3.00 


15.3 




2.80 


3.32 


3.45 


2.61 


3.89 


ii 


4.17 


18.3 


1.88 


3.97 


4.43 


4.85 


3.63 


5.47 


iii 


1.95 


10.0 




1.77 


2.22 


2.24 


1.68 


2.53 


iv 


2.69 


13.3 




2.48 


2.97 


3.09 


2.34 


3.89 


V 




3.78 














vi 


2.80 


14.6 




2.58 


3.09 


3.21 


2.41 


3.89 


vii 




3.93 















of the Higgs doublet can realize gauge coupling unification more naturally. Therefore, in 
the models v and vii it is not so easy to realize gauge coupling unification, though these 
models have the advantage that the FCNC constraint is weaker because ipi = t. Half 
integer charges can play the same role as Z2 parity or R-parity. For example, model i 
does not require R-parity, model vii does not require Z2 parity, and model vi requires only 
either of R-parity or Z2 parity. Models vi and vii have a value —1/2 for the charge of A, 
and therefore, because the unification scale is given by ~ A~ a , these models predict 
a more stable proton than the other models. Note that the model ii has the possibility 
that proton decay via dimension 5 operators dominates that via dimension 6 operators. 

4.1.2 E 6 scenario 

The model IE is the unique consistent model with a = — 1. Models I and I are charac- 
terized by a = —1/2 and the "effective charge" h e s = —17/4 and —4, respectively. These 
models of the E§ scenario require only either R-parity or Z2 parity, because half integer 
charges play the same role as these parities. Also, they automatically meet the condition 
for suppressing FCNC. 

4.2 Results 

For each model, we calculate the maximal value of the cutoff scale A max , with which gauge 
coupling unification is realized, setting the parameters y max , A and a s (M z ) as in Tables |] 
and |j. In the E§ scenario, we additionally adopt the constraint that the gauge coupling 
constants do not become too large. 

4.2.1 5O(10) scenario 

Roughly speaking, as h increases by 1, A max increase by a factor of approximately 1.5. 
On the other hand, for a single value of h (models i, iv and vi), the values of A max may 
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Table 4: A max in units of 10 16 GeV for some models of the Eq scenario. Basically, the 
parameters are set as y max = 2, A = 0.22 and aJ 1 (Mz) = 8.44, and aj(A^) < 1 is required, 
but in each case, one of these values is different from these, as indicated in the first row. 
The notation " — " means that there is no solution for gauge coupling unification with the 
conditions considered. Most cases with no solution result from the condition a^A^) < 1, 
which is much different from 5*0(10) cases.. 







v - 2 2 

y max ^ 


Z/max — 2 a 


A = 0.20 


A = 0.25 


a; 1 = 8.30 


a; 1 = 8.58 


1 loop 


I 


9.4 


79.0 


2.91 


4.57 


19.3 


10.2 


8.66 


11.0 


11 


11.7 


99.3 


3.63 


5.75 


21.1 


12.5 


10.7 


11.6 


m 




16.4 






6.10 






6.68 



differ from each other by only about 10%, and A max does not depend strongly on the 
parameters A and a s , as seen in Table |. If we set y max = 2, A max is not far from the 
naively expected value A ~ 2 x 10 16 GeV. Therefore, the prediction of the proton lifetime 
is at most a factor of ~ 20 larger than the naively estimated value 3 x 10 33 years for the 
models with a = —1, and 5 x 10 34 for the models with a = —1/2. Note that this factor 
of ~ 20 represents the largest value allowed by the ambiguity in the value of the 0(1) 
coefficients. This is an unlikely situation, and therefore we believe that the actual value 
of this factor is smaller. 

Unfortunately, the value A max is strongly dependent on y max , and unless we fix y max , 
we cannot precisely predict the lifetime of the proton. For this reason, we analyse more 
precisely the dependence on y max in Fig. [I]. It is found that the most probable A becomes 
smaller than the naively expected value A^ ~2 x 10 16 GeV, and this difference becomes 
larger as h decreases. Furthermore, it is obvious that gauge coupling unification is more 
difficult to realize for smaller h. Therefore, if we take account of the present limit on the 
proton lifetime provided by experiments, the models v and vii seem to be unrealistic. 

Finally, comparing the one- and two- loop results, it is seen that the two-loop effect in 
fact slightly worsens the situation, causing A max to decrease by a factor of ~ 1.4, which 
corresponds to a prediction of the proton lifetime that is shorter by a factor of ~ 1/4. 

4.2.2 Eq scenario 

Since the Eq scenario has a larger Higgs sector, the models in this case have more su- 
perheavy fields than in the SO (10) scenario. This means that the Eq models have larger 
degrees of freedom with regard to the values of the 0(1) coefficients than do the 5O(10) 
models. Therefore, it is expected that gauge coupling unification will be easier to realize 
in the Eq models. Actually, as seen in Table |], the maximal values of the cutoff A max 
tend to be larger here than in the 5*0(10) cases. However, since the Eq models have more 
superheavy fields, the gauge couplings around the unification scale tend to become so 
strong that the perturbative analysis is not reliable. For this reason, we add an additional 
constraint at the unification scale A^, 

a x (A A ) < 1. (1) 
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l0g 2 2/max 




i5.5 i6 16.5 17 " logA(GeV) 

Figure 1: y max vs A min and A max . Here, the four curves corresponding to the four models, i, 
ii, iii and vii {h = —3, —2, —4, —6, respectively), are plotted together. The allowed region 
is the upper region of the curve for each model. Other parameters are set as A = 0.22 
and a'^Mz) = 8.44. 

Mainly because of this constraint, the results displayed in Table £| behave significantly 
differently. In order to satisfy this condition that ax be small, for many cases in the E 6 
scenario, most of the degrees of freedom concerning the 0(1) coefficients are lost. We can 
understand this from the A dependence elucidated in Table [|. Because the parameter A 
represents a unit of mass of the superheavy fields, ax around the unification scale depends 
significantly on this parameter. Actually, as A becomes larger, the gauge couplings at 
the unification scale become smaller, and therefore more degrees of freedom of the 0(1) 
coefficients remain unfixed and can be freely adjusted for the purpose of realizing gauge 
coupling unification. Through this effect, the increase of the maximal value of the cutoff 
A max m the E e scenario becomes much larger than in the S'O(IO) scenario. In particular, 
in model H, ax becomes larger than in models I and E, and hence more degrees of freedom 
are lost due to the requirement that ax be small. Therefore it is more difficult to realize 
gauge coupling unification. Actually, when y max = 2, only the case with A = 0.25 can 
realize gauge coupling unification. If we tighten the constraint and require ax to be 
smaller than unity at the cutoff scale A, even in models I and E, coupling unification is 
difficult to realize. 

5 Discussion and summary 

In the 50(10) scenario, because the allowed range of values of the parameters A and 
a s (Mz) is limited, in this range, the maximal value of the cutoff A max does not change 
greatly. Because A max is obtained by tuning the 0{1) coefficients to give the largest value, 
a reasonable value of A would be less than the maximal value. If we fix y max = 2, the 



10 



maximal value of the cutoff does not differ by too much from the naively expected value, 
A ~ 2 x 10 16 GeV, and therefore a realistic prediction for proton decay can be considered 
to be not greatly different from the naive prediction, which is, for a = — 1, 

^--V,S Xl 033(^^) 4 (^«) 2 y e arS , (1 , 

and for a = —0.5, 

Unfortunately this prediction is strongly dependent on the unknown parameter for the 
(9(1) coefficients, y max . If y max = 4, the upper bound of the prediction may be beyond 
the scope of future experiments. 

In the Eq scenario, on the other hand, because the coupling constants tend to be 
too large to allow a reliable perturbative analysis, we have added an artificial constraint 
to preserve the validity of our perturbative analysis. Though the E e models have more 
degrees of freedom among the 0(1) coefficients than do the 5*0(10) models, most of 
them are lost to the requirement of suppressing the gauge couplings. If we had some 
mechanism to suppress the gauge coupling constants other than the degrees of freedom of 
the 0(1) coefficients, then the Eq models could much more easily realize gauge coupling 
unification. Actually, typical values of the maximal cutoff A max in the Eq cases are larger 
than in the 5*0 (10) cases. It may be more natural in the E 6 scenario to assume that the 
gauge couplings are in the non-perturbative region. Though the perturbative analysis is 
not reliable in such cases, we may expect that the actual result is not so different from the 
perturbative prediction. This follows from the argument that if gauge coupling unification 
in MSSM is not accidental, it may be natural to realize a situation similar to that in the 
perturbative region, which can explain gauge coupling unification in MSSM. 

Note that there are other effects that can change the running of the gauge couplings 
that are not taken into account in the analysis presented in this paper. One effect is 
from SUSY breaking parameters. Generically, the SUSY breaking scale is not one scale, 
and therefore there are some effects contributing to the condition for gauge coupling 
unification. Here, we roughly estimate the contribution in the case of model i of 50(10) 
GUT. If we set the masses of the colored superpartners to ITeV and the masses of colorless 
superpartners to 100 GeV, then A max decreases by a factor of ~ 0.4. It would seem that, 
this is probably an overestimation. For example, if we restrict the gaugino masses to 
satisfy the GUT relations 

M 1 = M 1= M 1 

0:3 Q.2 OL\ 

the factor by which A max is decreased becomes 0.5-0.6. The contribution of the SUSY 
breaking parameters tends to decrease A max . Therefore, the estimations given in this paper 
still give conservative limits on the proton lifetime. Another effect is from the lack of our 
knowledge about the values of the 0(1) coefficients. We have used 1 as the central value 
for the 0(1) coefficients, but this is only an assumption. The result necessarily depends 
on this central value, especially in E & case, because the mass spectrum of superheavy 
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fields depends on this value. We have not used the ambiguities of the gauge symmetry 
breaking scale due to the 0(1) coefficients. This effect can change not only the running of 
the gauge couplings directly but also the spectrum of superheavy fields if (CC^j > \~( C+C K 
From the non-renormalizable term XXCC, by developing a non- vanishing VEV, the field 
X can acquire a mass that is larger than expected. Though this effect exists only for 
fields lighter than A c+C , it may suppress the gauge couplings at the unification scale. 

Finally, we emphasize that even if the upper bound for the prediction of the proton 
lifetime is beyond the scope of future experiments, a more likely prediction of our analysis 
is near the naive prediction, which is not very far from the present experimental bound. 
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A Recipe 

As mentioned in §j|, in our scenario, the mass spectrum of every heavy Higgs is determined 
within the range allowed by the ambiguity of the 0(1) coefficients. Therefore we can easily 
calculate the second order (3 function of the effective theories appropriate for each scale 



(see Appendix |B|). If we use the DR scheme, [21] the naive step function approximation 
is good for connecting each gauge coupling constant of the neighboring effective theory, 
including the case in which the symmetries of these effective theories differ. 

We adopt SU(3) C x SU(2) L x U(1) Y for /i < A"^A, SU(3) C X SU(2) L x SU(2) R x 
U(1) B -l for A~^A < fx < \~ a A, SO(10) for A~ a A < /i < A in the 50(10) scenario and 
\~ a A < jj, < \~^~ A in the E 6 scenario, and E$ for A ~ A < /i < A in the E e scenario 
as the symmetries of the effective theories. Here, we approximate the masses of massive 
gauge bosons as being the same as the symmetry switching scale. 

Regarding the Yukawa coupling effect, we consider only that of the interaction directly 
related to the left- and right-handed top quark in terms of the relevant symmetry, i.e. 
U*Q 3 H U for SU(3) C x SU(2) L x U(1) Y , Q R3 Q L3 H D for SU(3) C x SU(2) L x SU(2) R x 
U(1) B -l, tf 3 (16)tf 3 (16)#(10) for 50(10), and ^ 3 (27)# 3 (27)$(27) for E 6 . For the scale 
dependence of the Yukawa coupling, we use one loop RGEs, which are sufficient for the 
purpose of investigating the flow of the gauge coupling constants. 

We start from the central values of the MS gauge coupling constants at \i = M z given 
in Ref . |p6| , including la ambiguity for a s (M z ). We set M 4 MS = 165 GeV, tan/3 = 5, 



and v = 174 GeV, which correspond to y t = 0.967. From Mz to the SUSY breaking 
scale, we use the RGEs of the standard model, which contain three family fermions and 
one Higgs doublet. Since we do not specify the SUSY breaking mechanism, we fix the 
SUSY breaking scale at 1 TeV and adopt a naive step function approximation, except 
in the transformation from the MS scheme into DR|27|], as the SUSY breaking threshold 
effect. At the SU(2) R x U{1)b-l breaking threshold, we must transform {ay, «2, a 3 } into 
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{«b-l, cur, 0L2-, a^}, in contrast to the one loop analysis. We set a]^_ L = ^(ay 1 — a^ 1 ) at 
/i = A~( C+C )/ 2 A, and iteratively adjust cxr so that it is equal to ct2 at the 5*0(10) breaking 
scale A~ a A. 

We investigate whether gauge coupling unification can be realized as explained above, 
by using the ambiguity of the 0(1) coefficients, for some values of the parameters A, 
2/max, A, and a s (Mz). In practice, we adjust all the independent masses of the superheavy 
Higgs by factors of y^ l l x -y mauX , instead of adjusting the 0(1) coefficients. Then, if all the 
differences of the values of a^ 1 can be made smaller than 0.05, we regard gauge coupling 
unification to be possible. 



B Renormalization group equations 

We use two-loop RGEs for the gauge coupling constants and one-loop RGEs for the top 
Yukawa coupling constant: 



0t 
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-b-+ ^ br+ M , 



16vr 2 1 (16vr 2 ) 2 v (16vr 



(1) 
(2) 



167T 2 167T 2 

ae constants. 

Considering only the top Yukawa interaction, as explained in Appendix |A], aj, C and 
Ci are determined as follows. |8| For SU(3) C X SU(2) L x U{l) Y (non-SUSY) (z = 1, 2, 3), 



Here, foj, 6^, a,, C and Cj are some constants. 



/ 17/10 

3/2 
V 2 



C = 9/2, C, 



17/20 
9/4 



(3) 



for 5C/(3) C x SU(2) L x [7(l) y (z = 1,2,3), 



26/5 
6 
4 



c = e, a 



13/15 
3 

16/3 



(4) 



for 5f/(3) c x SU(2) L x S77(2) fl x C/(1) B _ L (i — B — L, R, 2, 3), 



( 2 \ 

12 

12 
V 8 J 



c = 7, a 



/ 1/3 \ 
3 
3 

V 16/3 / 



(5) 



for 5O(10) (z = 5O(10)), 



and for i? fi (z =_E fi 



o< = 64, C = 40, C; = 63, 



(h = 180, C = 60, Ci = 104. 



(6) 
(7) 
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The constants bi and by are determined by the matter content of the effective theory. 
For example, for the standard model (i = 1, 2, 3), which contains three families of fermions 
and one Higgs doublet, 



( 41/10 \ 






kj = 







199/50 27/10 44/5 
9/10 35/6 12 
11/10 9/2 -26 



and for the MSSM (i = 1,2,3), which contains three family chiral superfields and two 
Higgs doublet chiral superfields, 



199/25 27/5 88/5 
9/5 25 24 
11/5 9 14 

The correction due to each representational chiral superfield is as follows. 
For SU(3) C x SU(2) L x U(1) Y (i = 1,2,3), 




A6,- 




1/75 3/5 16/15 
1/5 21 16 
2/15 6 68/3 

128/75 128/15 


16/15 34/3 

72/25 



8/75 32/15 


4/15 34/3 

/ 9/25 9/5 
3/5 7 

V o oo 










/ 


Ab l = 




, Ab i:j = 




12/ 


V 



21 16 



for Q + Q, 



for U c + U c 



for E c + E c 



for D c + D c 



for L + L, 



for G, 



for W, 



for X + X. 



(9) 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 

(17) 



For SU (3)c x SU (2) L x SU (2) rxU (1) b-Li first we define the abbreviation of each repre- 
sentation {SU{3) C , SU(2) L , SU(2) R ) u(1)b _ l as Q L (3, 2, l)i, Q R (3, 1, 2)_i, L L (1, 2, !)_!, 
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L R (1, 1, 2)i, H T (3, 1, 1)_2, [7(3, 1, 1)4, X Q (3, 2, 2)_ 2 and their conjugates, and H D (1, 2, 2) 

G(8, l,l)o, W L {1, 3, l)o, V fl (l, 1, 3)o and iV(l, 1, 1) . I The spinor 16, vector 10 and 
adjoint 45 of 50(10) are decomposed as 



o- 



16 -> Q L + Q R + L L + L R , 

10 -> H T + H T + H D , 

45 -> G + W^ + W^ + iV + Xr 



Xo + 17 + 17. 



Then, the corrections due to these representational chiral superfields are (i 
L, R, 2, 3) 



A6 ? ; 



/ 1/2 ^ 


/ 




3 


, Ab i:j = 


v 2 y 


\ 


/ 1/2 ^ 


I 


3 





v 2 y 







1/2 



1/2 








3/2 
21 



6 




21 
6 

! 





8/3 \ 



16 
68/3 / 

8/3 \ 
16 




Ah; 



V J 



Ah 



Ah 



/3/2\ 

1 


/ 3 / 2 \ 

1 







V 1 / 



A6,-,- 



68/3 / 

/ 9/4 9/2 ^ 



3/2 7 

V / 

/ 9/4 9/2 \ 



Ah 



3/2 





7 







0/ 



A6. 



/ 2/3 16/3 \ 













Ah 



Ah 



1 
1 

V j 

f°\ 




V 3 / 



Ah; 



Ah 



V 2/3 34/3 / 

/ ^ 

7 3 

3 7 

V / 

/ \ 





V 54 y 



>J 



for Ql + Ql, 



for Q R + Q R , 



for L L + L L , 



for L R + L R , 



for Ht + Ht, 



for H D , 



for G, 



(18) 
(19) 
(20) 

B - 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



3 Here we write U(V)b-l charges with the usual normalization, but in the context of SO(10) symmetry, 
they are multiplied by a/3/8. 
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Ah 



Ah 





2 

V / 

(°\ 

2 



V ) 



Ah 



( 4 \ 




V i J 





/ 


o 





\ 


Ahj = 

















24 






V o 








/ 




(° 








\ 


Ahj = 





24 




















I o 








/ 


/ 32/3 





64/3 


\ 



for W L , 



for Wr, 



Ab 



Ah 



( 4 \ 

6 

6 

V 4 y 



Ahj 






V 8/3 

/ 8/3 12 

4 42 

4 18 

V 8/3 12 






34/3 J 

12 64/3 \ 



for U + U, 



18 32 
42 32 
12 136/3 J 



forXn+Xr 



(2f 



(29) 



(30) 



(31) 
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